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Typeset using REVTgX 



Relativistic hadronic models suggest that nucleons propagating in nuclear matter develop 
large, attractive scalar and repulsive vector self-energies 0. Evidence supporting this picture 
in finite-density quantum chromodynamics (QCD) has emerged from some recent QCD 
sum rule calculations In particular, a rough connection is established between the 

nucleon scalar self-energy M* and the scalar condensate (qq) p in medium, and between the 
nucleon vector self-energy S v and the vector condensate (q'q) p - An alternative formulation 
of the sum rules ||, however, apparently gives qualitatively different results, in which the 
vector self-energy takes the opposite sign, so that the vector condensate induces attraction 
rather than repulsion. In this note, we resolve this ambiguity by properly accounting for 
spectral asymmetries in the finite-density sum rules for the nucleon, showing that the two 
formulations are each consistent with relativistic phenomenology. 

The QCD sum rules for a nucleon in nuclear matter focus on a correlator of interpolating 
fields for the nucleon: 

n N (g) = i J Jxe** ($ Q \T{ n (x)lj(0)}\*o) , (1) 

where |$o) is the nuclear matter ground state, which is characterized by four- velocity 
and baryon density p N in the rest frame. The interpolating field r](x) is taken here and in 
Refs. HH to be that advocated by Ioffe 0. (See Ref. |§ for a treatment of more general 
interpolating fields and for explicit expressions.) 

Because two four-vectors are available, the correlator can be decomposed into invariant 
functions in two ways: 1 

IT N (g) = ili q (q 2 } q ■ u) + U s (q 2 ,q ■ u) 

+ ill u (q 2 ,q-u) (2) 

= ffn 9 (g 2 ,g- «) + n s (g 2 ,g- «) 



1 Rclativistic covariance, time reversal, and parity imply that there are only three independent 
invariant functions 
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+ iiU u (q 2 ,q ■ u) , 



(3) 



where g M = g M — (q ■ u)u^. The two sets of invariant functions are trivially related: Tl q = Il 9 , 
Il s = Il s , and IL U = (q ■ u)U q + U u . These two alternative decompositions will lead to the 
two sets of sum rules we consider. 

The analytic structure of the correlator I1 N , and consequently the invariant functions 
in Eqs. (§) and is revealed by a standard Lehmann representation in energy lu, at 
fixed three-momentum q ||,|]|. [We work for convenience in the rest frame of the matter, 
where m m = (1,0) and g M = (u;,q).] The invariant functions are found to have singularities 
only on the real axis. In the upper half plane, IT N (u;,q) is equal to the retarded correlator 
and in the lower half plane to the advanced correlator. We will be concerned only with 
the discontinuities across the real axis (that is, the spectral functions), and therefore the 
infinitesimals that differentiate retarded, advanced, and time-ordered correlators for real uj 
are not relevant here. 

We exploit the analytic structure of the correlator by considering integrals over con- 
tours running above and below the real axis, and then closing in the upper and lower half 
planes, respectively (see Ref. @]). By approximating the correlator in the different regions 
of integration and applying Cauchy's theorem, we can derive a general class of sum rules, 
which manifest the duality between the physical hadronic spectrum and the spectral function 
calculated in a QCD expansion 



Here W(uS) is a smooth (analytic) weighting function and the spectral densities p ph and p th 
are proportional to the discontinuities of the invariant functions across the real axis. (These 
sum rules can also be derived by expanding dispersion relations for retarded and advanced 
correlators with external frequency u' in the limit to' — > zoo.) The phenomenological spectral 
density p ph models the low-energy physical spectrum, while the theoretical spectral density 
p th is approximated using an operator product expansion (OPE) applied to I1 N . The QCD 
sum rule approach assumes that, with suitable choices for W and the effective continuum 
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thresholds ooo and — cJ , each integral can be reliably calculated and meaningful results 
extracted. 

The thresholds ojq and — oJq act as effective boundaries beyond which the physical spec- 
trum p ph , when moderately smoothed, is identical to that of the leading contributions to p th . 
At zero density, the discrete space-time symmetries imply that p(— cu,q) = — p(+u;,q), and 
so we are led to take u7o = u . That is, charge conjugation implies that the free-space spec- 
tral function corresponding to positive-energy states (like the nucleon) is equal in magnitude 
to that of the corresponding negative-energy states (like the antinucleon). Since the relative 
sign of the positive- and negative-energy spectral functions is negative, only a weighting 
function odd in uj yields a nonzero result in Eq. (^). In this case, we can convert the integral 
to one over uj 2 and then over s = u 2 — q 2 . The end result features only manifestly covariant 
integrals, as one would expect. 

Useful choices for the weighting function at zero density are W(uj) = ue~ w2 ^ M2 , which 
generates the conventional Borel sum rules fTTH, or polynomials in uj, which generate the 
so-called finite energy sum rules (FESR) ]TD[ . (Actually, u times a monomial in s is used in 
the latter case, since p = p(s) in vacuum.) 

At finite density, ordinary charge conjugation symmetry is broken by the nonzero baryon 
number of the medium. This is clear from physical considerations: the propagation of a 
nucleon in ordinary nuclear matter is quite different from the propagation of an antinucleon. 
The consequence for the spectral function p ph (co>, q) is an asymmetry with respect to positive 
and negative u. Asymmetries also appear naturally in p th (uj, q) by generalizing the OPE to 
finite density [see below]. As a result, we are also compelled to adopt an asymmetric duality 
interval so that uo ^ uJo; this is the feature missed in Ref. @. 

Following Ref. we saturate the phenomenological integral in Eq. (^) with a quasipar- 
ticle pole ansatz: 

This implies 
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pf = ft 
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2E* 



X 2 5(lu - E q ) - \ 2 5(lu - E 



pf = Pf 
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M* 

2E* N 



t-2 , 



la 



q 

7T 



2E* q 
2E* q 



\ 2 5{u) - E g ) - \ 5(lu - E. 

\ 2 5{UJ - Eg) - \ 2 5{U - Eg) 

X 2 5(uj - Eg) + X 2 5(uj - Eg 



(6) 

(7) 
(8) 
(9) 



where E* = {q 2 + M* 2 ) 1 ^ 2 , E q = S v + E* and E q = £ 



E* < 0. The notation is that 



Pi (pi) corresponds to the discontinuity of IL; (Hi), with i = {q, s,u}. We have allowed for 
different residues A 2 and A for the positive and negative-energy poles, respectively. Different 
residues would naturally arise, for example, from the energy dependence of the self-energies. 
The form of the ansatz in Eq. (|5]) is otherwise constrained by Lorentz covariance. While a 
sharp quasinucleon state represents an extreme simplification of the actual spectrum, it is 
not unrealistic when smeared over energy scales of several hundred MeV, as we do in the 
sum rules. The physics motivation for the ansatz is discussed further in Ref. |Q. 

The operator product expansion for I1 N is developed in Refs. |4],[12]|| . We keep only 
the terms from the leading diagrams here (corresponding to those considered in the first 
paper of Ref. 0) to illustrate our point. For a quantitative analysis, one would include at 
least the contributions through terms involving four-quark condensates. The expansion is 
(suppressing parts that do not contribute to p th ) 

1 



n N (g) = i 



(<? 2 ) 2 in(-^ 



64tt 4 

+ :r^<? • uln(-q 2 )(q^q) 



3tt 2 



1 



^q 2 ln(-g 2 ) (qq) p + i^Q 2 ln(-<? 2 ) (q^q) 



P ! 



(10) 



which implies (specializing to the rest frame again) 



p q h = pf = n[9(co - \ci\ 



-u - q 



r 1 



{<ff 



1 



(11) 



x 



64vr 4 



3tt 2 



Ps 



th 




i h = -7r[%-|q|)-0(-a;-|q|)] 
,r[0(a,-|q|)-0(-a;-|q|)] 

\_3n z 



(12) 



,th 



(13) 




7r[fl(o; - |q|) - fl(-o; - |q|)] 



(14) 



We now have all the ingredients needed to construct QCD sum rules using Eq. (|J). 

The phenomenological expectation from Relativistic Brueckner Hartree-Fock (RBHF) 
and mean-field calculations, as well as from Dirac phenomenology, is that at nuclear matter 
density M* « 0.65M N and S v « 0.35M N , and the sum is relatively constant with density 



from higher-order condensates @||.) In these rules, M* is correlated with (qq) p , and S v > 
with (q ] q) p . 

The authors of two recent preprints, however, obtain very different qualitative results. 
In Ref. [|], the decomposition in Eq. @ was (implicitly) adopted, and applied with q = 0. 
The apparent advantage is that a direct projection onto the positive-energy quasinucleon is 
possible. That is, |(1 + 7o) projects onto the positive-energy quasinucleon and |(1 — 70) 
projects onto the negative-energy quasinucleon. From these projected sum rules, (q^q) p 
was found in Ref. || to be associated with attraction rather than repulsion, implying that 
S v < 0. The net result was a very large energy shift of the quasinucleon pole. Here we point 
out the missing elements in these calculations, which drastically affect the conclusions. 

As noted above, the essential point is that the asymmetry in the phenomenological 
spectral density must necessarily be reflected in an asymmetry in the duality interval; that 
is, UJq 7^ Uq. This point and its semi- quantitative implications are made clear by considering 
FESR's, using Eqs. (liTD-ffTi]) with monomials in 10 as the weighting functions W(u). 



|TJ. This expectation is supported by a simple version of the sum rules that is discussed in 
Ref. Q]. (A more complete sum rule is, however, not yet conclusive because of uncertainties 
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We start with the first decomposition, Eq. (|2|). Given a quasiparticle ansatz for the 
nucleon Eq. ([5]), this is a natural decomposition for isolating the scalar self-energy M* (by 
considering IE/ILJ and the vector self-energy S v (by considering H u /H q ). We apply Eq. (f|) 
to each of the three functions p q , p s , and p u , first with W{uS) = 1 and then with W{uS) = to. 
We can derive further sum rules, of course, but these are sufficient to make our point. For 
simplicity and for a clear comparison to Ref . || , we take q = from here on, so that q 2 — > uo 2 
in Eqs. ©-(13) ■ We obtain 

1 :(A 2 - A 2 ) = - or*) - ^q)M + tf) , (15) 



u 



2M* X ' 320vr 4V u uy 6tt 2 

^[M*(A 2 + A 2 ) + £ V (A 2 - A 2 )] = + &o) - 9^ + ?»o - &S) , (16) 

l^-* 2 ) = -^mM-ini), (I?) 

1[M*(A 2 + A 2 ) + £ V (A 2 - A 2 )] = ~^- 2 {qq)M + D$) , (18) 

^(A 2 -A 2 )4aK 3 -^), (19) 



2M ;[K(A 2 + A 2 ) + S V (A 2 - A 2 )] = ^(^K + ^j . (20) 

One should be cautious about quantitative results extracted from these sum rules, but the 
qualitative features should persist in more sophisticated treatments. 

The terms in the operator product expansion (OPE) odd in u vanish at zero density. 
But at finite density, they imply, through duality, an asymmetry in the physical spectral 
functions, as manifested in Eq. fll5D . There are many more such terms as we extend the 
OPE to higher dimension. 

If we expand in density, then from Eq. (^) we see that (recall that A 2 = A and uj = uJ 
at p N = from charge conjugation symmetry) 

u -ul ~ 0(pjs) , (21) 
A 2 -A 2 ~0(p N ). (22) 

If we drop terms of higher order in p N , Eqs. (0), fll8| ), and (|20|) become 
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^(A 2 + A 2 ) = ^iK 6 +^), (23) 
iM*(A 2 + A 2 ) = - I ^(gg) P K 4 + ^) , (24) 
^ V (A 2 + A 2 ) = ^(gt g)p (^ + ^). (25) 



We can divide Eqs. (E§ and (031) by (E3) to find 



24tt 2 647r 2 
M n = ~t~2\^p and Sv = Tir(? ?)p ' ( 26 ) 

where (a; 2 ) ~ i(a>o + 57q) to this or der. Plugging in a typical continuum threshold we 
find roughly the magnitudes, density dependence, and scalar-vector cancellations expected 
for the self-energies from relativistic phenomenology. Note that S v is unambiguously pos- 
itive. The scale of u — oJq is set by Eq. fllBD, which in conjunction with Eq. fllTf ) and the 
results for S v and M* in Eq. fl26|), implies that lj — ZJ rs S v . 

Now we consider the second decomposition, Eq. fl3|), and extract the analogous sum rules. 
The rules from IT^ and IT S are the same as above. The only different sum rules come from 

flu- 

u: ^ A2 + ^V^K 6 +^)-^K 3 -^), (27) 

|(A 2 + A 2 ) + ^(A 2 - A 2 ) = -^K 4 + O + ^(-o 7 " *5) • (28) 

The underlying problem with Ref . is made clear by Eq. ( p8"D . If uJq = u is assumed [and 
therefore A = A 2 from Eq. (fPTDl, one concludes that S v < 0. 

Adding the sum rules for f[ s = U s and TL U (that is, making the |(1 + 70) projection), we 
obtain: 

A2 = 385?<"» + - IS^ -=>!>- -=* . (29) 

A*(E, + K) = -^f + - + + J^ll-o 7 " 4) ■ (30) 

Dividing these equations gives a sum rule for the positive-energy pole: 



E q (q = 0) = S v + M* 



-iMp + ±(q*q) P )H + uj a ) + - gg) 



(31) 



This is the FESR analog of the Borel sum-rule equations derived in Ref. ||. Again, if 
ZJq = uo is assumed, (q^q) p and the change with density of (qq) p both apparently reduce E q 
from its zero density value of M N . However, once uJq ^ ujq is allowed, several new terms 
contribute. Numerical estimates using ujq — ujq = S v from Eq. (|26|), condensate values from 
Ref. [|J, and ojq ~ 1.6-1.7 GeV show that E q (q = 0) is essentially constant with density 
rather than decreasing precipitously, as found in Ref. || .0 

This quantitative result, which involves various cancellations, is very sensitive to details 
of the calculation: the choice of weighting function, the continuum threshold, the number 
of terms kept in the OPE. In contrast, sum rules derived from the first decomposition of 
I1 N , which lead to separate expressions for the scalar and vector self-energies rather than 
for the sum, are more robust. In Refs. the decomposition of Eq. (0) was considered, 

with an asymmetric weighting function used to suppress contributions from the negative- 
energy quasinucleon in favor of the positive-energy quasinucleon. The possibility u ^ uJ 
was not considered in these references; however, explicit calculations show that making 
this generalization is numerically unimportant. Thus there is a qualitative difference in the 
sensitivity of the two correlator decompositions [Eqs. (0),(|J)] to the (uncertain) details of 
the continuum contribution. 

In summary, we have reexamined the results of Ref. ||, which implied that certain QCD 
sum rules predict that the vector self-energy of a nucleon in nuclear matter is attractive, 
rather than repulsive as implied by relativistic phenomenology and other sum rules j|. 
By properly accounting for asymmetries in the spectral functions between positive- and 
negative-energy states, we have shown that there is no inconsistency, and that all formula- 
tions lead to a repulsive vector interaction. 

2 These rather high values for ujq are needed in the FESR at p N = to obtain a sensible result for 
M N . In addition, we note that the sum rules imply that ujq +£Jo is roughly constant with density. 
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